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Statement  of  Problems  Studied 


This  report  summarizes  the  research  carried  out  on  the  development  of  a  numerical 
approach  for  simulating  radiative  emission  under  conditions  of  thermochemical  nonequilibrium 
that  are  produced  in  high  altitude,  hypersonic  flows. 

Summary  of  Most  Important  Results 

This  report  summarizes  the  research  carried  out  on  the  development  of  a  numerical 
approach  for  simulating  radiative  emission  under  conditions  of  thermochemical  nonequilibrium 
that  are  produced  in  high  altitude,  hypersonic  flows.  We  have  undertaken  the  development  of  a 
hybrid  computer  simulation  method  for  these  flows  that  uses  a  continuum  approach  based  on 
traditional  Computational  Fluid  Dynamics  (CFD)  solving  the  Navier-Stokes  equations,  and  a 
kinetic  approach  based  on  kinetic  theory  using  the  direct  simulation  Monte  Carlo  method 
(DSMC).  Prior  to  this  study,  no  such  hybrid  code  existed  in  the  literature  for  hypersonic  flows. 
The  report  included  as  an  Appendix  describes  the  first  version  of  our  hybrid  approach  that  uses 
the  Information  Preservation  (IP)  method  as  the  communication  link  between  the  CFD  and 
DSMC  methods.  The  hybrid  code  is  applied  to  a  Mach  4  flow  over  a  wedge.  Comparisons 
between  CFD,  DSMC,  and  hybrid  methods  reveal  strengths  and  weaknesses  of  the  present 
approach. 
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Assessment  of  a  Hybrid  Method  for  Hypersonic  Flows 

Wen-Lan  Wang*,  Quantum  Sun*  and  Iain  D.  Boyd* 


*  Department  of  Aerospace  Engineering ,  University  of  Michigan,  Ann  Arbor,  Ml  48109 


Abstract  A  preliminary  hybrid  particle-continuum  computational  framework  for  simulating  hypersonic  interacting  flows 
is  proposed*  The  framework  consists  of  the  direct  simulation  Monte  Carlo-Information  Preservation  (DSMC-IP)  method 
coupled  with  a  Navier-Stokes  solver.  Since  the  DSMC-IP  method  prorides  the  macroscopic  information  in  each  time 
step,  determination  of  the  continuum  fluxes  across  the  interface  between  the  particle  and  continuum  domains  becomes 
straightforward.  A  hypersonic  flow  over  a  two-dimensional  wedge  is  considered  as  an  example  and  compared  with  pure 
particle  calculations.  The  results  show  that  this  preliminary  hybrid  framework  is  promising  but  several  issues  are  yet  to  be 
resolved. 


INTRODUCTION 

In  recent  years,  biconic  configurations  and  hollow  cylinder/flares  were  extensively  studied  with  different  numerical 
approaches  [1, 2,  3, 4]  for  a  set  of  hypersonic  free-stream  condition.  In  comparisons  with  experimental  data  [5],  it  was 
evident  that  the  computational  fluid  dynamics  (CFD)  techniques  based  on  the  Navier-Stokes  (NS)  equations  have  better 
performance  than  the  direct  simulation  Monte  Carlo  (DSMC)  method  [6]  in  terms  of  the  capture  of  flow  structures  and 
the  prediction  of  the  separation  region  size.  However,  It  was  shown  that  the  CFD  methods  constantly  over  predict  the 
heat  transfer  on  the  fore  body  ahead  of  where  separation  occurs.  Furthermore,  it  was  demonstrated  [7]  that  adjacent 
to  the  body  wall,  especially  near  the  cone  tip,  there  is  a  strong  thermal  nonequilibrium  effect,  that  is  inconsistent 
with  the  basic  assumptions  of  the  NS  equations.  In  addition,  Roy  et  al.  [8]  showed  that  DSMC  and  CFD  do  agree 
for  conditions  that  are  more  rarefied  than  the  experiments.  Therefore,  the  objective  of  the  present  work  is  to  develop 
a  combined  DSMC-CFD  computational  method  for  the  physically  accurate  and  numerically  efficient  analysis  of  this 
class  of  hypersonic  flows. 

Two  primary  Issues  must  be  taken  into  account  in  developing  a  combination  of  DSMC  and  CFD  methods.  We  first 
need  to  know  when  to  switch  between  the  methods.  Since  it  is  well  known  that  the  NS  equations  are  not  valid  under 
rarefied  conditions,  it  is  general  to  use  a  continuum  breakdown  parameter  as  the  criterion.  For  the  hypersonic  flows 
mentioned  above,  this  issue  has  been  investigated  in  a  previous  study  [9]  and  where  it  was  concluded  that  a  proposed 
parameter 

Knmax  =  max  (Kn^,  Kn^?  Kn^)  (1) 

with  a  threshold  value  of  0.05  can  best  predict  the  regions  where  the  Navier-Stokes  equations  fail.  The  Knudsen  number 
in  Eq.  1  is  expressed  as  Kn^  =  g  |  VQ|,  with  Q  being  any  flow  property  In  general.  Since  continuum  breakdown  is  often 
related  to  the  transport  phenomena  of  viscosity  and  heat  transfer,  we  only  consider  the  flow  properties  of  density(D), 
temperature(T)  and  magnitude  of  velocity(V). 

Another  issue  regarding  the  development  of  a  hybrid  method  is  the  approach  for  information  exchange  at  the 
Interface  between  the  DSMC  and  CFD  domains.  At  the  interface,  macroscopic  flow  properties  must  be  provided  to  the 
CFD  method  to  evaluate  the  net  fluxes  and  to  the  DSMC  method  to  initialize  the  particles  entering  from  the  continuum 
domain  Into  the  rarefaction  domain.  Several  approaches  have  been  considered,  such  as  the  Marshak  condition  [10], 
the  kinetic  flux- vector  splitting  (KFVS)  scheme  [11,  12]  and  the  adaptive  mesh  and  algorithm  refinement  (AMAR) 
embedding  a  particle  method  [13].  Unfortunately,  the  development  of  a  robust,  multi-dimensional  scheme  that  is 
capable  of  handling  nonequilibrium,  hypersonic  compressed  flows  has  not  yet  been  accomplished. 

The  primary  difficulty  of  the  second  issue  results  from  the  fact  that  the  DSMC  approach  always  involves  very  strong 
statistical  scatter  unless  the  number  of  samples  Is  large  enough.  Recently,  the  information  preservation  method  (IP) 
was  proposed  to  reduce  the  statistical  scatter  in  low-speed,  constant  density  flow  systems  [14].  Since  then,  the  method 


has  been  generalized  to  allow  density  and  temperature  variations  and  shown  to  be  effective  in  solving  microfluidic 
and  micro-electro-mechanical  system  (MEMS)  problems  [15,  16,  17].  One  advantage  of  the  IP  method  is  that  the 
macroscopic  values  of  the  flow  field  are  known  at  any  time,  since  the  information  in  cells  is  updated  for  each  time  step. 
This  eases  the  complexity  of  coupling  the  particle  method  with  the  continuum  solver.  A  hybrid  approach  that  combines 
the  IP  method  and  a  NS  solver  has  shown  great  progress  in  solving  micro-scale  gas  flows  [18].  This  motivates  us  in 
the  present  study  to  explore  the  IP  method  for  hypersonic,  interacting  flows. 

The  layout  of  the  paper  is  as  follows,  A  brief  description  of  the  continuum  approach  will  be  presented  first  in  the  next 
section,  followed  by  an  introduction  of  the  DSMC-IP  method,  A  detailed  explanation  on  how  to  combine  these  two 
methods  together  is  provided  in  the  section  of  Domain  Coupling,  A  hypersonic  flow  over  a  two-dimensional  wedge 
is  considered  to  assess  the  new  hybrid  technique.  In  the  last  section,  conclusions  and  suggestions  for  future  work  are 
provided. 


NUMERICAL  SCHEMES 
Continuum  Approach 

The  Navier-Stokes  equations  in  the  continuum  domain  are  solved  numerically  using  an  explicit  Gauss-Seidel  line 
relaxation  method  and  second-order  accurate,  modified  Steger- Warming  flux  vector  splitting  [19],  The  viscosity  p  is 
modeled  with  the  power  law  and  the  thermal  conductivity  rc  is  determined  from  the  Prandtl  number 

^l  =  ^ref(r/rref)e>5  ^  =  cpVl/k 

where  cp  is  the  specific  heat  at  constant  pressure.  A  slip-boundary  model  proposed  by  Gokgen  [20]  is  implemented. 


DSMC-IP  Approach 


The  information  preservation  method  was  first  developed  by  Fan  and  Shen  [14]  to  overcome  the  statistical  scatter 
problem  in  DSMC  simulations,  especially  for  systems  in  which  the  flow  speed  is  much  smaller  than  the  molecular 
speed.  In  addition  to  the  ordinary  thermal  velocity  that  is  utilized  to  compute  the  particle  trajectory,  each  simulation 
particle  in  the  DSMC-IP  method  also  possesses  macroscopic  preserved  information  such  as  velocity  vector  and 
temperature.  The  DSMC-IP  method  has  achieved  great  success  for  solving  micro-scale  gas  flows  [see  15, 16, 17]. 

In  the  most  recent  work  by  Sun  and  Boyd  [21],  an  additional  temperature  term  is  introduced  in  order  to  solve 
the  contradiction  between  the  real  flux  and  the  DSMC-IP  representation  of  the  translational  energy  flux  across  a  cell 
interface.  From  gas  kinetic  theory,  the  average  translational  energy  of  a  molecule  at  equilibrium  temperature  T  is  3 kT /2 
(k  is  the  Boltzmann  constant),  whereas  the  average  translational  energy  carried  by  a  molecule  across  an  interface  is 
2kT,  The  extra  energy  must  be  correctly  modeled  or  the  energy  of  the  whole  system  will  not  be  conserved.  Therefore, 
the  simulation  particles  in  the  IP  method  have  an  additional  temperature  Ta. 

In  each  time  step  of  the  DSMC-IP  method,  simulation  particles  are  first  moved  and  collided  in  the  usual  way  as 
in  the  standard  DSMC  method.  The  preserved  velocity  in  the  ri  direction  and  temperature  of  simulation  particles  are 
updated  by  solving 
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where  p  is  the  pressure,  p  is  the  mass  density,  £  is  the  number  of  internal  degrees  of  freedom  of  molecules,  and 
the  subscript  c  denotes  the  macroscopic  information  for  the  computational  cells.  After  the  preserved  information  of 
simulation  particles  is  updated,  the  preserved  information  for  cells  is  updated  by  taking  the  arithmetic  average  over  the 
information  of  all  Np  particles  in  the  cell. 
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The  density  is  updated  by  solving  the  continuity  equation 
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The  ideal  gas  law,  p  =  pRT,  is  assumed.  A  detailed  description  and  implementation  of  the  DSMC-IP  method  can  be 
found  in  Ref.  [21]. 

Note  that  an  adequate  numerical  scheme  must  be  employed  to  solve  Eq.  6  due  to  the  presence  of  shock  waves  in 
supersonic  flows.  Since  the  continuity  equation  also  appears  in  the  NS  equations,  we  use  the  same  technique  described 
in  the  last  subsection  for  solving  it. 

The  current  DSMC-IP  code  is  based  on  a  parallel  optimized  DSMC  code  named  MONACO  [22].  A  sub-cell  scheme 
is  implemented  for  selection  of  collision  pairs  where  the  number  of  sub-cells  is  scaled  by  the  local  mean  free  path. 


DOMAIN  COUPLING 

To  implement  the  coupling  between  the  particle  method  and  the  NS  solver,  buffer  and  reservoir  DSMC-IP  cells  are 
introduced  in  the  continuum  domain  adjacent  to  the  domain  interface,  as  depicted  in  Fig.  1 .  A  similar  concept  of 
reservoir  cells  was  first  proposed  in  Ref.  [23]. 

The  buffer  DSMC-IP  cells  work  as  an  extension  of  the  particle  domain.  Simulation  particles  that  end  their  movement 
phase  within  the  pure  particle  domain  or  in  the  buffer  cells  are  retained.  Those  that  leave  these  two  regions  are  removed. 
For  each  time  step,  all  simulation  particles  in  the  reservoir  cells  are  first  deleted  and  then  re-generated  based  on  the 
cell-centered  values.  The  number  of  new  particles  is  evaluated  from  the  cell  density  value  and  the  particle  velocities 
and  temperature  are  initialized  to  the  Chapman-Enskog  distribution  [24]  based  on  the  corresponding  cell  values.  The 
newly  generated  particles  are  randomly  distributed  within  the  reservoir  cells.  In  this  study,  one  layer  of  buffer  cells  and 
five  layers  of  reservoir  cells  are  employed. 

In  the  continuum  domain,  the  NS  solver  determines  the  interface  continuum  fluxes  by  using  the  NS  variables  and 
DSMC-IP  cell  macroscopic  information.  Since  the  macroscopic  information  in  the  DSMC-IP  cells  is  known  in  each 
time  step,  the  DSMC-IP  cells  adjacent  to  the  domain  interface  are  treated  as  the  ghost  cells  that  provide  the  boundary 
conditions  for  flux  computations. 


NUMERICAL  EXAMPLE 

This  section  describes  a  Mach  4  numerical  experiment  to  assess  the  hybrid  technique.  A  geometry  of  a  2D  wedge  with 
a  25°  half  angle  is  chosen.  The  symmetric  line  of  the  wedge  is  aligned  with  the  free-stream.  Therefore,  only  the  upper 
half  of  the  wedge  is  considered.  Since  the  wake  region  behind  the  wedge  is  not  of  interest  in  this  investigation,  we 
assume  that  the  wedge  is  infinitely  long  but  only  the  first  5  cm  from  the  leading  edge  is  considered.  A  structured  grid, 
300  cells  along  the  body  by  200  cells  normal  to  the  body,  used  in  all  our  computations  is  shown  in  Fig.  2. 

The  fluid  is  pure  nitrogen  and  the  free-stream  conditions  are:  £/*»  =  1111.1  m/s.Too  —  1S5.6  K  and  pfX3  =  6.545  x  10~4 
kg/m3.  The  mean-free-path  in  the  free-stream  is  about  10-4  m.  An  isothermal  wall  at  a  temperature  of  293.3  K  is 
assumed.  These  specific  flow  conditions  are  the  same  as  those  in  CUBRC  Run  28  [25]  except  for  a  lower  flow  speed. 

In  the  CFD  calculation,  jU^  =  1 .656  x  105  N*s/m2,  I ^  =  273  K  and  (D  =  0.74.  The  Prandtl  number  is  considered 
as  a  constant  of  0,72. 

A  pure  DSMC  steady  state  solution  is  obtained  with  the  use  of  more  than  1.8  million  simulation  particles  at  the  end 
of  the  computation.  The  reference  time  step  in  the  pure  DSMC  calculation  is  5  nsec.  300,000  time  steps  of  computation 
are  performed  and  the  last  50,000  time  steps  are  sampled  to  obtain  the  results. 

In  the  hybrid  simulation,  a  steady  state  solution  from  the  CFD  method  is  first  obtained.  Using  the  steady  state 
solution  and  Knmax,  one  can  determine  the  locations  of  the  interfaces  between  the  continuum  and  the  particle  domains. 
In  the  particle  region,  cell  values  are  set  to  the  CFD  steady  state  results  and  simulation  particles  are  initialized  to  the 
Chapman-Enskog  distribution. 


FIGURE  1.  Interface  cell  types. 


FIGURE  3,  Particle  domain  in  hybrid  simulation:  (a)  initial  particle  domain  and  (b)  particle  domain  after  400,000  iterations. 


Although  it  has  been  concluded  in  a  previous  study  [9]  that  Kn max  =  0.05  is  a  good  criterion  for  separating  the  two 
domains,  a  more  conservative  criterion  of  Knmax  =  0.02  is  actually  employed  in  this  investigation.  The  initial  particle 
domain  and  the  domain  after  400,000  hybrid  simulation  iterations  are  illustrated  in  the  shaded  areas  of  Figs.  3.  About 
40.2%  of  the  total  number  of  cells  are  in  the  particle  domain  at  the  beginning.  This  value  increases  gradually  with  time 
and  reaches  about  44.8%  in  the  end.  Each  cell  in  the  particle  domain  has  about  70  simulation  particles  on  average.  The 
reference  time  step  in  the  hybrid  computation  is  1  nsec.  The  hybrid  method  results  presented  In  this  paper  are  obtained 
by  sampling  the  last  20,000  time  steps. 


Flow  Fields  and  Surface  Properties 

Comparisons  of  the  density  and  translational  temperature  contours  obtained  with  pure  DSMC  and  the  hybrid 
methods  are  made  in  Figs.  4(a)  and  4(b),  respectively.  The  hybrid  method  results  are  in  good  agreement  with  the 
pure  DSMC  results,  in  general.  Additionally,  the  hybrid  method  provides  much  smoother  results. 

In  Fig.  5,  comparisons  of  flow  properties  along  the  wedge  surface  are  made.  Generally  speaking,  the  flow  properties 
predicted  by  the  hybrid  method  along  the  wall  are  in  excellent  agreement  with  the  DSMC  solutions.  It  is  interesting 
to  note  that  the  heat  transfer  rate  (St =qw/  jpooU^)  predicted  by  the  hybrid  method  slightly  shifts  away  from  the  initial 
CFD  result  that  is  almost  identical  to  the  DSMC  result.  A  similar  behavior  is  also  observed  for  the  skin  friction 
coefficient  in  Fig.  5(b).  In  contrast  to  heat  transfer  rate  and  skin  friction,  pressure  along  the  wedge  surface  computed 
with  the  hybrid  method  moves  towards  the  DSMC  result,  as  shown  in  Fig.  5(c). 


FIGURE  4,  Comparison  of  DSMC  and  hybrid  solutions  of  (a)  mass  density  and  (b)  translational  temperature. 


FIGURE  5.  Comparison  of  surface  properties  with  different  numerical  approaches:  (a)  heat  transfer  rate,  (b)  skin  friction 
coefficient  and  (c)  surface  pressure. 


Detailed  Comparisons 

In  Fig.  6,  the  pure  DSMC,  CFD  and  hybrid  solutions  are  compared  in  detail  along  a  line  normal  to  the  wedge  at 
x  ~  0.5  mm.  Since  this  location  is  very  close  to  the  wedge  leading  edge,  where  the  flow  is  strongly  nonequilibrium,  the 
CFD  and  pure  DSMC  solutions  are  therefore  different  to  a  large  extent.  It  is  evident  that  in  the  region  very  close  to  the 
wall  (Sn  <0.15  mm),  the  density  and  velocity  obtained  with  the  hybrid  approach  are  in  outstanding  agreement  with 
the  pure  DSMC  results.  Overall,  they  both  move  from  the  initial  CFD  conditions  towards  the  DSMC  solutions  and 
are  qualitatively  consistent  with  the  DSMC  solutions.  However,  the  translational  temperature  predicted  by  the  hybrid 
approach  is  very  different  to  the  DSMC  results.  The  reason  for  this  poor  result  may  come  from  the  translational  energy 
model  and  will  be  discussed  later.  The  horizontal  dotted  line  in  Fig.  6  indicates  the  interface  between  the  particle 
region  (lower  portion)  and  continuum  region  (upper  portion). 

A  similar  comparison  between  the  pure  DSMC,  CFD  and  hybrid  solutions  for  density,  velocity  and  translational 
temperature  at  x  —  1  mm  is  displayed  in  Fig.  7.  It  is  clear  that  the  hybrid  technique  works  fairly  well  in  computing  the 
density  and  velocity  but  poorly  in  modeling  the  translational  temperature. 

Further  downstream  at  x  =  5  mm  where  the  nonequilibrium  effect  starts  to  weaken,  the  Navier-Stokes  equations 
coupled  with  a  slip-boundary  model  perform  very  well.  As  shown  in  Fig.  8,  the  CFD  and  DSMC  solutions  approach 
one  another,  especially  in  the  region  adjacent  to  the  wall. 

The  whole  computational  domain  at  this  location  is  divided  into  four  regions  along  the  Sn  direction:  particle  (I: 
0  <  Sn  <  0.8),  continuum  (II:  0.8  <  Sn  <  1.86),  particle  (III:  1.86  <  Sn  <  2.6)  and  continuum  (IV:  5n  >  2.6).  It  is 
shown  in  Fig.  8(a)  that  the  density  and  translational  temperature  predicted  by  the  hybrid  approach  are  in  fair  agreement 
with  the  pure  DSMC  solutions,  except  near  the  interface  between  regions  II  and  III.  The  velocity  profile  in  Fig.  8(b) 
exhibits  unexpected  behavior  at  the  interface  between  regions  I  and  II.  The  reasons  for  that  may  be  due  to  the  size 


FIGURE  6.  (a)  Temperature  and  density  and  (b)  velocity  profiles  along  the  line  normal  to  the  wedge  surface  at  x  =  0.5  mm. 


FIGURE  7.  (a)  Temperature  and  density  and  (b)  velocity  profiles  along  the  line  normal  to  the  wedge  surface  at  x  =  1  mm. 


of  region  III  being  too  small  or  the  numerical  scheme  employed  to  solve  Eq.  6  in  the  particle  regions,  and  will  be 
discussed  later.  Similar  comparisons  at  x  =  2  cm  are  plotted  in  Fig.  9. 


Discussion 

There  are  three  major  issues  regarding  the  hybrid  method  results.  First,  it  is  clear  that  the  breakdown  criterion 
Knmax  =  0,02  does  not  capture  the  shock  front  very  accurately,  causing  the  particle  region  near  the  shock  to  be  too 
small  and  in  turn  causing  the  hybrid  technique  to  have  difficulties  to  accurately  predict  the  shock  structures.  As  a  result, 
it  is  difficult  for  the  continuum  method  after  the  shock  to  obtain  the  correct  solution  because  the  DSMC-IP  method 
provides  a  questionable  boundary  condition. 

The  second  issue  occurs  in  the  region  very  close  to  the  wedge  leading  edge  where  the  translational  temperature 
is  not  correctly  predicted  by  the  DSMC-IP  method.  As  mentioned  in  the  Numerical  Schemes  section,  the  DSMC- 
IP  method  has  an  additional  temperature  term  for  balancing  the  energy  transported  by  a  molecule  across  an  interface. 
However,  this  extra  energy  estimated  in  the  current  DSMC-IP  method  is  based  on  the  equilibrium  assumption  in  which 
the  velocity  distribution  function  is  Maxwellian.  Because  of  the  strong  nonequilibrium  near  the  leading  edge,  the  extra 
energy  modeled  in  the  DSMC-IP  technique  is  no  longer  valid.  Consequently,  the  translational  temperature  is  poorly 
predicted  by  the  DSMC-IP  approach  in  the  particle  region. 

The  hybrid  velocity  profiles  in  Figs,  8  and  9  is  the  third  concern.  The  reason  for  the  unexpected  behavior  may  be 
attributed  to  the  first  issue,  or  it  may  be  a  result  of  other  factors,  such  as  the  way  we  solve  the  continuity  equation 
Eq.  6  in  the  particle  region.  In  this  study,  we  adopt  the  same  numerical  scheme  as  used  in  the  CFD  code  to  solve  Eq.  6. 


FIGURE  8.  (a)  Temperature  and  density  and  (b)  velocity  profiles  along  the  line  normal  to  the  wedge  surface  at  x  =  5  mm. 


FIGURE  9.  (a)  Temperature  and  density  and  (b)  velocity  profiles  along  the  line  normal  to  the  wedge  surface  at  x  =  2  cm. 


On  the  cell  interfaces,  the  way  that  the  specific  numerical  scheme  in  CFD  calculates  the  mass  flux  largely  depends 
on  the  characteristic  equations  which  are  obtained  from  the  Euler  equations.  By  definition,  neither  the  Euler  nor  the 
Navier-Stokes  equations  are  valid  in  the  continuum-breakdown  region.  As  a  result,  using  any  techniques  derived  from 
them  may  induce  errors. 


CONCLUDING  REMARKS  AND  FUTURE  WORK 

We  have  described  a  preliminary  hybrid  particle-continuum  method  for  computing  hypersonic,  interacting  flows.  In 
the  particle  domain,  the  DSMC-IP  technique  was  employed  while  in  the  continuum  domain,  a  finite  volume,  second- 
order  accurate,  Steger- Warming  flux  vector  splitting  NS  solver  was  utilized.  Since  the  macroscopic  values  of  the  flow 
field  in  the  DSMC-IP  approach  are  updated  in  cells  for  each  time  step  advanced,  coupling  these  two  domain  becomes 
straightforward.  Information  exchange  on  the  interface  between  the  particle  and  continuum  domains  was  continuously 
carried  out  for  each  time  step. 

A  numerical  experiment  of  a  hypersonic  flow  over  a  wedge  was  conducted,  with  the  use  of  the  continuum-breakdown 
criterion  Knmax  =  0.02.  The  hybrid  method  was  first  initialized  with  a  steady  state  CFD  solution  and  then  marched 
forward.  It  is  concluded  that  the  new  hybrid  method  is  able  to  simulate  the  hypersonic  flow  very  well,  generally 
speaking.  The  detailed  comparisons  showed  that  the  density  and  velocity  computed  with  the  hybrid  method  matched 
with  the  exact  (pure  DSMC)  solutions  over  the  entire  flow  field  faithfully,  but  translational  temperature  was  poorly 
modeled  near  the  leading  edge. 


The  poor  performance  of  the  hybrid  method  in  predicting  the  translational  temperature  around  the  tip  of  the  wedge  is 
primarily  caused  by  the  equilibrium  approximation  in  the  DSMC-IP  technique  to  estimate  the  extra  translational  energy 
carried  by  a  molecule  moving  from  one  cell  to  the  other.  A  variant  of  this  equilibrium  model  is  under  investigation. 

An  adequate  numerical  scheme  for  solving  the  continuity  equation  in  the  particle  region  is  necessary  as  the  Euler 
and  Navi er- Stokes  equations  fail  in  the  continuum-breakdown  region  and  consequently  any  numerical  schemes  derived 
from  these  equations  may  fail  in  this  region. 

The  numerical  experiment  also  showed  that  the  size  of  the  particle  region  around  the  shock  wave  is  too  small  using 
the  current  breakdown  criterion.  Some  portion  of  the  shock  front  was  not  included  in  the  particle  domain.  This  aspect 
may  lead  to  the  unexpected  behavior  observed  in  the  velocity  profile.  A  more  conservative  criterion  would  help  to 
solve  this  problem  but  would  also  reduce  the  numerical  efficiency. 
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